Quark electric dipole moment induced by magnetic field 
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We show numerically that quarks develop an electric dipole moment in the direction of a suf- 
ficiently intense magnetic field due to local fluctuations of topological charge. This anomalous 
CP-odd effect is a spin analogue of the chiral magnetic effect in QCD. 
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I. INTRODUCTION 

Fluctuations of topological charge is an important fea- 
ture of the QCD vacuum which is reflected in the mass 
spectrum of hadrons, in the breaking of chiral symmetry, 
and other properties. Recently it was suggested that the 
topological fluctuations may be directly observed in hot 
QCD matter in the presence of a very intense (hadron- 
scale) external magnetic field via the so called "chiral 
magnetic effect" (CME) The fluctuations of the topo- 
logical charge lead to a local imbalance between the left- 
and right-handed light quarks. In the presence of an 
intense magnetic field the quarks move along the field, 
and the chiral imbalance creates a net electric current 
along the direction of the magnetic field. As a result, 
this CP-odd effect should lead to a separation of the 
electric charges along the field, which potentially may be 
observed in experiment. 

The strong magnetic fields of the order of QCD scale 
may be created in noncentral heavy-ion collisions due 
to relative motion of electrically charged ions in an ini- 
tial state and the products of the collision in a final 
state [l| Q . The strong fields are perpendicular to the 
collision plane, so that the CME should be realized on an 
event- by- event basis in a form of an asymmetry of elec- 
trically charged particles emitted below and above the 
reaction plane. There is a recent experimental evidence 
reported by the STAR collaboration that this scenario 
works at RHIC [3] . Recently, the CME was also found in 
numerical lattice simulations 4] (a brief review of lattice 
QCD simulations in strong magnetic fields can be found 
in Ref. [1). 

In this paper we show that the topological fluctuations 
in the background of the intense magnetic field lead to 
appearance of a quark electric dipole moment. This effect 
is a spin analogue of the CME. Indeed, the CME induces 
an electric dipole moment due to the spatial (nonlocal) 
separation of the electric charges of quarks while the dis- 
cussed spin effect induces the electric dipole moment lo- 
cally. 

In Sec. fll] we discuss general features of this spin effect. 
Sec. Illll is devoted to a numerical demonstration of the ef- 
fect in a simple case of an instanton-like configuration. In 



Sec. lIVI the quenched lattice simulations are used to show 
that the quark electric dipole moment is induced both in 
high- and low-temperature phases. Our conclusions are 
given in the last section, and the appendix contains tech- 
nical details of our simulations. 



II. GENERAL ARGUMENTS 

Quarks are the electrically charged spin-1/2 particles 
which carry certain magnetic dipole moments. A strong 
enough external magnetic field aligns quark spins and 
leads to appearance of a net magnetization. A quantita- 
tive measure of the magnetization is given by the expec- 
tation value [6| 



where 



S Q /3 = Yi^ al13 ~ ll3la 



(1) 



(2) 



is the relativistic spin operator, F a p = d a ap — dpa a is the 
strength tensor of the external electromagnetic field, and 
a a is the Abelian gauge field. For the sake of simplicity 
we consider one quark flavor, and omit the flavor and 
spinor indices in Eq. (JT|). 

The structure of the right hand side of Eq. ([TJ is obvi- 
ous: the appearance of the electromagnetic field strength 
tensor F a p is dictated by a covariance argument, while 
the appearance of the chiral condensate (3 , ^ r ) is the nat- 
ural consequence of the dimensionality of the left and 
right parts of Eq. (JT|). Moreover, the left hand side of ([T]) 
violates the chiral symmetry similarly to the chiral con- 
densate. The coefficient of proportionality x{F) i s called 
the chiral magnetic susceptibility @ . 

There is another property which was used for the 
parametrization of the right hand side of Eq. (JXJ) : this 
equation is formulated in the CP-symmetric vacuum. If 
the vacuum is not CP-symmetric, then Eq. (JXJ) should be 
rewritten in a more general form: 

(*£a/J*>fl = (**) X(F) gF aP + x(F) qF a0 , (3) 
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where 



F, 



a(3 — c^otfiiivF^u 



(4) 



is the dual field strength tensor. In Eq. ^ we used the 
subscript (• • • )g to stress that the spin expectation value 
is evaluated in a CP-odd environment, which can be in- 
duced, say, by a nonzero 6>-angle. We omit the subscript 
9 in the right hand side of Eq. © . The new coefficient 
X is the electric susceptibility of the vacuum which in 
general may be nonzero if 9 =/= 0. 

One can rewrite Eq. ([3]) by introducing the dimension- 
less magnetic and electric dipole moments, respectively 

Ma/3 = X( F ) qF afj , e a/3 = x(F) qF a/3 , (5) 

so that Eq. (J3J is simplified: 

(*S Q(3 *) e = + e aP ) (tf *) . (6) 

Here the tensors (i a p and e Q/ g represent two orthogonal 
contributions, 



Pap £a/3 



0. 



(7) 



Note that in presence of a nonzero 9- angle the ten- 
sor structure of the magnetization ([3J is most general. 
The nonzero (9-angle leads to appearance of the pseu- 
doscalar condensate (^75 the pseudotensor conden- 
sate (^>j 5 t a G^^), where G£„ is the gluon field strength 
tensor and t a are the generators of the gauge group, and 
other condensates. The pseudoscalar condensate may 
give separate contributions to the chiral susceptibilities 
X and x m Eq. (|3j, while it cannot mix the two terms 
in ([3]) because of their mutual orthogonality. As for the 
pseudotensor condensate, it may be split - in terms of the 
Lorentz structures - into two orthogonal terms that are 
proportional to the electromagnetic field strength ten- 
sor F^ v and its dual (QJ, respectively. Therefore the 
pseudotensor condensate may also contribute to the both 
terms of Eq. <J3]) leaving the general form of the magne- 
tization (PJ unchanged. 

It is convenient to rewrite the magnetic and electric 
dipole moments ([5]) via the expectation values ([3]) in the 
vector form, respectively: 



(^ 75 S i0 ^) = 1 

m) - 2 £ijk 
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(8) 
(9) 



where we used the relation 75 E Q ^ = £ Qj a, and assumed 
that the external electric field is absent, Foi = 0. The 
interpretation of Eq. © as the electric dipole moment 
is obvious because this vector quantity is defined by the 
electric component of the spin operator, EjQ. Here spatial 
tensor indices are i, j = 1, 2, 3 and the Euclidean time is 
labeled by the index 0. 



In the real CP-invariant vacuum the electric dipole 
moment © should be identically equal to zero. However, 
the fluctuations of the electric magnetic dipole moment 
may still be very strong. The real vacuum may be sub- 
divided into topologically nontrivial CP-odd domains in 
which the quarks may have possess an anomalous elec- 
tric dipole moment along the direction of the magnetic 
field, e || B. In these domains the CP-odd electric dipole 
moment may be as strong as the magnetic dipole mo- 
ment. Below we check this conjecture using numerical 
simulations. 

In order to characterize quantitatively the appearance 
of the quark electric dipole moment, we consider the fol- 
lowing quantities: the local magnetic dipole moment 
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the local electric dipole moment 

trf (x) = $ (x) EioV" (x) , 
and the local chirality 

P5 (x) = $ (x) 7 5-0 (x) = p L (x) - p R (x) 



(10) 



(11) 



(12) 



which is the operator of the difference of the densities of 
left- and right- handed quarks. 



III. INSTANTON 

Instanton is a topologically nontrivial solution to the 
classical equations of Yang-Mills theory. Topological 
charge of the instanton is nonzero, and therefore it can be 
considered as a simplest gauge filed configuration which 
may lead to the anomalous electric dipole moment. 
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FIG. 1: The local density of the magnetic dipole moment of 
quarks (|10|l in the instanton background exposed to a strong 
magnetic field. The upper (blue) and lower (green) surfaces 
represent, respectively, the longitudinal, a^ 1 (x), and a trans- 
verse, o"i f (x), components of the dipole moment in a 34- plane. 

While a spectrum of the Dirac operator in the field 
of a single instanton is well-known [7], the presence of 
a uniform magnetic field makes the problem analytically 
intractable. Therefore we constructed an instantonlike 
smooth configuration of the lattice gauge field with unit 
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FIG. 2: The same as in Fig. [JJ but for the local density of the 
electric dipole moment 



topological charge and numerically calculated fermionic 
propagators in such a configuration (details of our nu- 
merical simulations are given in the appendix). Next, we 
added the uniform magnetic field in the third direction, 
Bi = BSiz t° this sample configuration. We evaluated the 
local density af 1 of the magnetic dipole moment (|10[) in 
the background of this configuration. 

Note that any gauge field configuration - including the 
instanton - has, in general, both zero and near-zero Dirac 
eigenmodes. These classes of the Dirac modes have dif- 
ferent physical meaning. For example, the number of 
the zero modes is related to the topological charge of the 
gauge field configuration according to the Atiyah-Singer 
theorem. Consequently, the vacuum expectation value of 
the number of zero modes is proportional to the suscep- 
tibility of the topological charge. 

On the contrary, the Banks-Casher relation [8| states 
that the vacuum expectation value of the density of the 
near-zero modes is proportional to a different quantity, 
the chiral condensate. Therefore, the near-zero modes 
carry information about the chiral properties of the sys- 
tem, while the exact zero mode, in general, does not. In 
Ref. [10( we have analytically related the chiral magne- 
tization to the spin structure of the near-zero modes in 
the external magnetic field (we discuss this relation in 
the Appendix, and refer an interested reader to Ref. [To[ 
for more details). Thus, only near-zero modes were taken 
into account in our calculations of the chiral magnetiza- 
tion. 

We would like to make a remark that in various instan- 
ton gas and liquid models the individual zero modes of 
instantons and antiinstantons do play an important role 
in the chiral properties of the system because in this in- 
teracting system the zero-mode degeneracy is lifted out 
and the zero modes become eventually near-zero. In our 
illustrative example we consider the single instanton like 
configuration so that zero mode does not contribute to 
the chiral magnetization. 

The density of the longitudinal (i = 3) and the trans- 
verse (i = 1) components of the magnetic dipole moment 
in the 34-plane of the background instantonlike configu- 
ration are plotted in Fig. [TJ The distribution in the 12- 
plane is similar to the one plotted in Fig. [TJ As one can 
expect, the magnetic dipole moment is predominantly 



directed along the magnetic field: Icrf^x)! 3> l^rSC 1 )! 
so that a M \\ B with a good accuracy. Notice that the 
spatial distribution of the density of the magnetic dipole 
moment is uniform similarly to the magnetic field. 

The density of the electric dipole moment is shown in 
Fig. [3J Similarly to the magnetic moment, the electric 
dipole moment ([TTj) is directed along the magnetic field: 
I cr E (x) I 3> |cf 2 (a;)| so that a E || B with a good accuracy. 
The magnitude of the electric dipole moment is smaller 
and, in contrast to the magnetic dipole moment, its dis- 
tribution is peaked near the position of the instanton. 

Summarizing this section, we observed the genera- 
tion of the electric dipole moment in the background of 
the topologically nontrivial configuration of the smooth 
gauge fields. In the next section we study the same effect 
in a real ground state of non-Abelian gauge theory. 

IV. QUENCHED QCD 

A. Fluctuations of electric dipole moment 

The vacuum expectation value of the magnetic dipole 
moment of the quark (|10l) is given by (UJ: 

(§ct m *) = X {B) qB . (13) 

The chiral magnetic susceptibility x was calculated nu- 
merically in Ref. [l(| • 

The average of the electric dipole moment (fTTj) is zero, 

<*ct e *)=0, (14) 

since the density of topological charge changes its sign in 
different space-time domains, thus supporting the global 
CP-invariance of the vacuum. However, locally the elec- 
tric dipole moment a E may be nonzero, as we have seen 
in the previous section for the instanton case. 

In order to measure the strength of the local fluctu- 
ations of the dipole moments we study the connected 
expectation values 

<W) 2 >,« = (W-«» 2 ) BiT 

-<W-<«0)*> W . < 15 > 

where there is no summation over the indices i and £. 
In Eq. (|15p erf is the i-th component of the magnetic 
(£ = M) or electric (I — E) dipole moment, and (. . .) B T 
denotes the expectation value with respect to the ther- 
mal state at temperature T in the background magnetic 
field B. The subtraction of the expectation value at 
B = 0, T = removes ultraviolet divergences and yields 
physical results which are practically independent of the 
UV cutoff (this question was discussed in Refs. 0,[I(| us- 
ing the same set of the gauge field configurations). The 
subscript IR Eq. (|15l) reflects the fact that the subtrac- 
tion provides us with the nonperturbative infrared (IR) 
value. 
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FIG. 3: The fluctuations (|15[) of the longitudinal (i.e., directed 
along the magnetic field) components of the magnetic (|10p and 
electric dipole densities vs the strength of the magnetic 
field in the confinement (T = 0) and in the deconfinement 
(T = 1.12 T c ) phases. 



We calculated the fluctuations (fT5|) numerically in 
SU (2) lattice gauge theory with quenched massless chi- 
rally invariant quarks (details of simulations can be found 
in the appendix). In Fig.[3]we show the fluctuations (fT5)) 
of the longitudinal components of the electric dipole 
moment (0 = , Eq. (jTTJ)) both in the confinement 
phase, at T = 0, and in the deconfinement phase, at 
T = 1.12 T c . Here T c rj 310 MeV is the critical temper- 
ature of the confinement-deconfinement phase transition 
in the SU(2) gauge theory. For comparison, we also plot- 
ted the longitudinal component of the magnetic dipole 
moment [O — a^ 1 , Eq. (fTU|) ]. One can see that the aver- 
age squares of the magnetic moment and of the electric 
dipole moment are equal with a good precision. Both of 
them grow significantly with magnetic field. Taking into 
account the fact that the average electric dipole moment 
is zero, one can conclude that there are space-time do- 
mains within which the magnetic dipole moment and the 
electric dipole moment are either parallel or antiparallel, 
and are equally strong. It is interesting that the fluc- 
tuations of both magnetic and electric dipole moments 
are almost independent of the temperature even near the 
phase transition (the small mutual deviations in Fig. [3] 
are of the order of the error bars) . 

We also found that the fluctuations of the transverse 
components of the electric af 2 an d magnetic a^ 2 dipole 
moments are almost independent of the magnetic field. 
For these components the expectation value (fT5j) is com- 
patible with zero. 

Summarizing this subsection, we conclude that in the 
external magnetic field the quark develops the electric 
dipole moment, which is directed along the magnetic 
field. The fluctuations of the electric dipole moment are 
of the same order as the magnetic ones. 



B. Electric dipole moment and chirality 

In order to demonstrate that the electric dipole mo- 
ment is closely related with the local chirality we calcu- 
late the correlator of the electric dipole moment with the 
chiral density (fT2")l : 



( e,m\ 

[P5,^i ) 



I E.M\ 



(16) 



In this formula no summation over the index i is implied. 
We plot the correlator (fT6|) of the longitudinal electric 



dipole moment, cf' , in Fig. U both in the confinement 
phase and in the deconfinement phase. At zero tempera- 
ture the correlator grows quickly and we observe almost 
the full (100%) correlation of the quark electric dipole 
moment and the chirality even at weakest nonzero mag- 
netic field. The strength of the effect is somewhat smaller 
in the deconfinement phase at weak magnetic fields due to 
thermal fluctuations. However, at strong magnetic fields, 
qB ~ 1 GeV 2 , the correlation function (|16l) reaches the 
highest possible value, c (ps, erf') = 1. 
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FIG. 4: The correlation (|16[) of the longitudinal component 
of the electric dipole moment (llip of quark (EDM) with the 
chiral density (|12|l vs the strength of the magnetic field in 
the confinement and deconfinement phases. The horizontal 
dotted line indicates the 100% correlation. 

The transverse components af 2 of the electric dipole 
moment (fTTj) are not locally correlated with the chiral- 
ity (fT2|) . The same is true both for transverse a^ 2 an d 
longitudinal a^ 1 components of the magnetic dipole mo- 
ment of quark. Note that no correlation is observed in 
the absence of the external magnetic field. 

Summarizing, at zero temperature the electric dipole 
moment is strongly correlated with the chiral density in 
the presence of the external magnetic field. The ther- 
mal fluctuations reduce this correlation, which is restored 
again as the strength of the magnetic field increases. 
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V. CONCLUSION 

We demonstrated that the quark develops the anoma- 
lous electric dipole moment in the presence of a suffi- 
ciently strong external magnetic field. The anomalous 
moment is parallel to the direction of the magnetic field. 
The observed effect is a reflection of a topological struc- 
ture of the QCD vacuum: the induced electric dipole 
moment is strongly correlated with the chiral density of 
quarks which, in turn, is correlated with the topological 
charge density. The sign of the electric dipole moment 
of the quark is a fluctuating quantity so that the electric 
dipole moment of the quark is zero on average, and no 
global CP-violation occurs. However, the local fluctu- 
ations of the longitudinal (i.e., parallel to the magnetic 
field) component of the quark electric dipole moment are 
as strong as the fluctuations of the conventional magnetic 
dipole moment. 

The effect is observed at a single instanton gauge field 
configuration, as well as in a true vacuum state of a 
quenched non-Abelian gauge theory. The effect is strong 
both in the confinement and deconfinement phases. The 
thermal fluctuations decrease the correlation of the lon- 
gitudinal electric dipole moment with the chiral density 
at moderate magnetic fields. As the strength of the field 
increases, these two quantities become fully correlated. 



Appendix A: Details of simulations 

The setup of our numerical simulations is the same as 
the one used in Refs. 0,113, EH- We utilize lattice QCD 
with the simplest SU(2) gauge group because the gen- 
eration of the electric dipole moment originates in the 
chiral sector of QCD and therefore the number of colors 
is not crucial. In our simulations only valence quarks in- 
teract with the electromagnetic field. The effects of the 
virtual quarks on gluons are neglected because the in- 
clusion of dynamical (sea) quarks makes the simulations 
computationally difficult, while the essential features of 
the studied effect remains intact in the quenched limit. 

In Table|T]we present the parameters of our lattice sim- 
ulations: the lattice geometry, L^L t , the coupling con- 
stants /?, the lattice spacings a, the spatial lattice exten- 
sion L = L s a, and the minimal value of the magnetic 



fields \/qB n 



TABLE I: Parameters of simulations. 



T/T c L 3 L t P a, fm L s a, fm 



MeV 



14 14 3.281 0.103 1.44 343 
1.12 16 6 3.325 0.095 1.33 371 



for the gluon fields (see, e.g., Eq. (1) in 13]). The uniform 
magnetic field B in the direction \x — 3 is introduced into 
the Dirac operator by substituting su (2)-valued vector 
potential with u (2)-valued field: ->■ A% + , 
where C M = B/2 (a^c^i — xiS^)- Notice that both su(2) 
and u(l) are algebra-valued functions. 

This expression is valid in the infinite volume. In order 
to make it consistent with the periodic boundary condi- 
tions in the spatial directions, we have introduced an 
additional ^-dependent boundary twist for fermions on 
the lattice [lij . The finiteness of the volume leads to the 
quantization of the total magnetic flux, so that 



qB = 2tt k/L 2 



k G Z, 



(Al) 



where q = 1/3 e is the smallest (absolute value of) electric 
charge of the quark, and L is the length of the lattice in 
the spatial direction. 

In order to calculate the chiral expectation values |T|), 
(fT5]) and ([To]) we use the basis of the eigenmodes ip n of 
the Dirac operator T> = 7 M (<9 M — iA^), 



T>tp n = \ n ip v 



p(X) = {J2S{X-X n )), (A2) 



where A„ are the eigenvalues of the Dirac operator. 

The chiral condensate is calculated using the Banks- 
Casher formula M, 



lim lim 

A->0 V->oc 



MA) 
V 



(A3) 



where V is the total four- volume of Euclidean space-time, 
and p (A) is the density of the Dirac eigenvalues. In the 
quenched approximation the averaging is performed over 
the gauge fields A^ with the weight exp (— SVm [^D, 
where SVm [A^] is Yang-Mills action. In the Euclidean 
space the spinor conjugation is given by the complex con- 
jugation ip a = ipl l . 

In |10| we derived a magnetization analogue of the 
Banks-Casher formula (IA3I): 



which is used for evaluation of the expectation values con- 
taining the spin operator. The expectation values involv- 
ing four fermionic fields are evaluated using the following 
formula: 



(ipO^ipOuip) = Tr 



' Ox I Tr [ — — 2 



-Tr 



V + m 
1 



1 



D + m T> + m 



V + m 



2 ), 



In order to implement chirally symmetric massless 
fermions on the lattice, we use Neuberger's overlap Dirac 
operator [l2| ]. Ultraviolet lattice artifacts are reduced 
with the help of the tadpole-improved Symanzik action 



where 0\ and O2 are some spinor operators. The re- 
sult is then averaged over all configurations of the gauge 
fields. The Dirac propagator is evaluated by inverting 
the massive Dirac operator in the subspace spanned by 
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M Dirac eigenvectors which correspond to M nonzero 
Dirac eigenvalues with smallest absolute values: 



1 



V 



m 



= Yl 



k<M 



4>k (x) 4> k (y) 

Afe + m 



(A4) 



The value of M is limited by the numerical proce- 
dure used to hnd the eigensystem of the Dirac opera- 
tor (ARPACK in our case). We have used M = 10 and 
M = 12 (after the subtraction of ultraviolet divergences 
the expectation values are almost independent of M for 
M > 10, Ref. (3). 
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